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Effect of Spin Current on Uniform Ferromagnetism: Domain Nucleation 
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Large spin current applied to a uniform ferromagnet leads to a spin-wave instability as pointed 
out recently. In this paper, it is shown that such spin-wave instability is absent in a state containing 
a domain wall, which indicates that nucleation of magnetic domains occurs above a certain critical 
spin current. This scenario is supported also by an explicit energy comparison of the two states 
under spin current. 
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Magnetization dynamics driven by spin torque from a 
spin polarized current (spin current) has been studied ex- 
tensively after the theoretical prediction jlj, ^ that a spin 
current can be used to flip the magnetization in pillar (or 
spin valve) structures 0,SIESSI3- As a theoretical 
framework to describe such current-induced magnetiza- 
tion dynamics, Bazaliy, Jones and Zhang (BJZ) 0| de- 
rived a modified Landau-Lifshitz-Gilbert (LLG) equation 
for a fully-polarized ferromagnet (half metal) with a new 
term proportional to j ■ Vn. This new term represents 
a spin torque that a current of density j exerts on a lo- 
cal magnetization n. Later on, it was generalized to the 
partially polarized case with an interpretation of j as the 
spin current density, js 0, ^3 ■ In the Hamiltonian, the 
effect of spin torque appears in a form 
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where 9 and are polar angles which parameterize n. As 
seen from this form, spin current favors a magnetic con- 
figuration with spatial gradient, or more precisely, with 
finite Berry-phase curvature. It is thus expected that a 
large spin current destabilizes a uniform ferromagnetic 
state. This is indeed seen from the spin-wave energy 
illlllll, 
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where A — ^ J /K and n — K^/K, with J, K and K±^ 
being exchange constant, easy- and hard- axis anisotropy 
constants, respectively, of localized spins. The first term 
is the well-known spin-wave dispersion with anisotropy 
gap. The effect of spin current appears in the second 



term as the "Doppler shift" where i's(oc js) rep- 

resents drift velocity of electron spins. For sufficiently 
large v^, fi™' becomes negative for a range of k. This 
means that there exist states with negative excitation en- 
ergy, indicating the instability of the assumed uniformly- 
magnetized state (9(,|l^|llj|. The true ground state under 
spin current, however, remains to be identified. 

In this Letter, we point out that possible ground state 
is a state containing domain walls. Namely, the energy of 
a domain wall becomes lower than the uniform ferromag- 
netic state when a spin current exceeds a certain critical 
value . The behavior of nucleated domains depends 
much on the magnetic anisotropy parameters; Flowing 
domain wall is stable when the anisotropy is of uniaxial 
type. For the case of large hard-axis anisotropy, static 
domain wall is stable at the nucleation threshold, but 
when it starts to flow under larger current, even a state 
with a domain wall becomes unstable. The nature of the 
resulting state is still unknown. 

Our prediction of domain formation by spin current 
may be related to the very recent experimental obser- 
vations in metallic and semiconducting pillars and films 
[I1.1I1IIHIEIIE1I3, which suggest spatially inhomoge- 
neous magnetization reversal. Domain nucleation is also 
suggested by recent numerical simulation • 

We start by extending the formulation of BJZ to 
arbitrary degree of polarization, and derive the effective 
Lagrangian for magnetization which is slowly varying in 
space and time by assumption. With this effective La- 
grangian, we calculate spin-wave dispersion around a do- 
main wall solution in the presence of spin current. The 
spin-wave Doppler shift term now has the form, k{Vs—X), 
where X is the speed of the domain wall. From this re- 
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suit, we find that the spin-wave instability does not occur 
around a domain wall for any large spin current if the 
anisotropy is of uniaxial type {K > K±^). This suggests 
that the true ground state is a multi-domain state with 
domain walls. These domain walls turn out to be flowing 
with average speed equal to drift velocity Vs of electron 
spin determined by the spin current. The stability is then 
understood from Galilean invariance, since a domain wall 
moving with speed X = is equivalent to a domain wall 
at rest and in the absence of spin current. 

The situation is slightly different in the case of large 
hard-axis anisotropy, K± » K. Domain walls created 
by the spin current above j" are at rest if the current 
is below depinning threshold, j^opm oc K±X{> j"). At 
higher current, the walls start to flow, but this triggers 
another spin-wave instability. Thus multi-domain state 
collapses for js > into another new ground state 

which is still unknown We defer this new state to 
future studies. 

We consider a ferromagnet consisting of localized spins 
and conduction electrons. The spins are assumed to have 
an easy z-axis and a hard y-axis, and described by the 
Lagrangian 



Ls = 
Hs = 
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We have adopted a continuum description for localized 
spins, S{x,t) — Sn{x,t), with unit vector n{x,t) — 
(sin 9 cos 0, sin 9 sin (jj, cos 9) and the magnitude of spin, 
S. The J is the exchange constant, and a the lattice con- 
stant. The easy-axis (K) and hard- axis {K±) anisotropy 
constants generally incorporate the effect of demagne- 
tizing field. The exchange interaction between localized 
spins and conduction electrons is given by 



A / d^xn{x,t)-{c'{x,t)(Tc{x,t)), 
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where 2A and c (c^) are the energy splitting and anni- 
hilation (creation) operator of conduction electrons, re- 
spectively, and cr is a Pauli-matrix vector. The free- 
electron part is given by = J^ka '^kci^Cka with 
Efc = h^k^ /2m. We perform a local gauge transforma- 
tion [23, 1^ in electron spin space so that the quanti- 
zation axis is parallel to >S'(a;, t) at each point of space 
and time; c{x,t) — U{x,t)d{x,t), where a{x,t) is the 
two-component electron operator in the rotated frame, 
and U{x,t) = m{x,t) ■ cr is an SU(2) matrix with 
m(x,t) = (sin(6'/2)cos0,sin(6'/2)sin(/),cos(6'/2)). The 
electron part of the Hamiltonian is now given by Hc\ = 



teraction Hamiltonian by 
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Here, A^{q,t) = ^ d^xe"'-'^''^ A^{x,t) with A^(a;,t) = 
m X df^m represents SU(2) gauge field with space (^ = 
x,y, z) and time (/i — 0) components. 

For slowly varying magnetic configurations, the effec- 
tive Lagrangian can be derived by a perturbative expan- 
sion with respect to A^. This is the gradient expan- 
sion. In conformity with Ls, we retain terms up to sec- 
ond order in spatial gradient or in Ai (i = x,y,z), and 
first order in time derivative or in Aq. The spin-torque 
term, ■ , arises from the first-order contribu- 

tion {Hint). Here js = f J2ka ^^fka is the spin current 
density in the rotated frame, with a distribution func- 
tion /fcCT = (aL(i)afecr(i)) specifying the state carrying 
current. Other terms just give renormalizations of J and 
S ;20j . The effective Lagrangian is thus given by 



LcS ~ La, — Hi 
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The second term, given by Eq. J^l, with js is identical 
to that derived by BJZ ;9] for the half-metallic case. It is 
seen that spin current favors a finite Berry phase curva- 
ture along the current. The modified LLG equation can 
be obtained by taking variations of Leff with respect to 9 
and (/). It has the same form as in BJZ, and is applicable 
to general slowly-varying spin texture, and to the case of 
arbitrary degree of polarization. 

Let us now study the spin-wave excitation around a 
domain-wall solution under spin current by using the 
method of collective coordinate jl^ . The wall is assumed 
to be planar and moving in the direction perpendicular to 
the plane (chosen as x-direction) . It is convenient to use 
a complex field £,{x,t) = e"^*-^'*^ tan to represent 

magnetization. We decompose it into the domain-wall 
configuration and spin waves around it: 



£,{x,t) = ^dw(2:,t)exp 



2cosh{(x-X)/A} 



y^^ipk(,x - X)rik{t) 



(8) 



where ^dw = exp[— (x — X{t))/\ + i(f)o{t)\ represents a 
domain- wall configuration. The wall position X{t) and 
the angle (j)o{t) between wall spins and the easy plane 
are proper collective coordinates 0, H^, |3 EE H^- 
The spin- wave part is expanded with modes ^Pk{x) = 
-^=^={-ik\ + tfmh'j^)e''''' mill, whose amplitude 
rjk precisely corresponds to the Holstein-Primakoff bo- 
Sfccr ^kaO-k^cika with Efccr = h^k^ /2m — uA, and the in- son. (L is the length of the system in the z-direction.) 
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Substituting ((SJ into Q, we obtain the spin- wave La- 
grangian isw, up to the quadratic order in rj, as 



+_£(! _ 4sin2 (l,o){v-kVk + V*-kVl)j, (9) 

where N = 2XA/a^ is the number of spins in the waU {A 
is the cross-sectional area of the wall), fi^"' = {k^X^ + 
I + k/2) and Vc = ^^^^ ■ Details of the calculation will 
be presented in a separate paper. From Eq. @, we find 
the spin-wave dispersion around a current-driven domain 
wall: 

1/2 



-I- 2K^sin-'0oCOs20o| + k [vs - X j . {10) 

The most important difference from the case of uniform 
ferromagnet is that the drift velocity Vs — ^5 Js of spin 
current appears as a relative velocity, Vg — X, with respect 
to the moving wall. 

Equation H10|l shows that the spin-wave energy de- 
pends strongly on the wall dynamics (through X and (t>n) ■ 
From the equations of motion for the domain wall |2lj, 
A</io + C(X = and X — aX(f>o = Vc sin2(/io -f- Vs, where 
a <C 1 is a damping parameter, we have {1 + a'^)X = 
Vc sin 200 + Vs- Neglecting terms of ^ 0{a^), we have. 
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which depends on js implicitly through 0o • 

For K > K±, it is easy to see that fi^^ remains pos- 
itive for all values of k and js |2^ . Thus a domain- wall 
state is stable irrespective of the magnitude of the spin 
current in this case. This indicates that a uniform fer- 
romagnctism collapses into a multi-domain configuration 
for js > if, where jf = ^KX{1 + VTT^) is the crit- 
ical spin current density for the instability of a uniform 
ferromagnet _Q, IQ^ JlJ . Since the domain wall flows for 
js > jf"'^''' = t^KXk 113, which satisfies > p''^, it 
starts to flow as soon as it is created (Fig. 1(a)). The 
stability of moving domain wall state is natural from 
Galilean invariance since domain wall with velocity of 
Us is equivalent to the static domain wall without spin 
current. 

In the opposite case, K± > K , the spin wave shows no 
anomaly as long as js < j^''^'", where the wall remains 
static with constant 0o (= 



for jf < js < jf (jf < jf^'P''' is realized when 
K± > 8K), the static multi-domain state is realized. 
In contrast, as soon as the wall starts to move under 
larger spin current (jg > j^''''™), the spin wave becomes 
unstable (Eg. ((111) '). Thus another ground state would ap- 
pear. The wavelength of the unstable mode around the 
uniform magnetization, k = —{K{K + _ftr^)/ J^)^/^[llJ, 
suggests that this unknown ground state has a spatial 
structure with short length scale of ~ {KK±/ J^)^^^"^ = 
{K/K±y^^X < X. 

We have thus shown that the spin- wave instability un- 
der spin current is avoided by the existence of domain 
wall. This may indicate that domain nucleation occurs 
under spin current. We cannot, however, answer here 
how many domains are created. The above argument 
holds for each segment larger than the wall thickness. A, 
and thus the domain wall can be nucleated with spacing 
of ^ A. For a correct estimate of the spacing, however, we 
need to consider the dipolar interaction among domains. 
We will not pursue this point in this letter. 

We here present a supporting argument for the above 
scenario of domain formation by evaluating explicitly 
the total energy of a domain wall in the presence of 
spin current. We consider only the case K <C K± and 
j" < js js'^^™, where the wall can remain static and 
thus the energy comparison has physical meaning. A 
static domain wall configuration is given by 6{x,t) = 
9s{x) = 2tan~-^ exp(— a;/A'), (j){x,t) = 4's{x), where A' = 
i)o includes the effect of contraction of do- 
'1/ dx(j)s{dxds) siuOs). Note that 



A/ vl + Ksii? 
main wall 28] (0o 
(j)s{x) is not spatially constant when there is a spin cur- 
rent. For large hard-axis anisotropy, K± '3> -^jjs, sub- 
stitution of the static solution 9s and 0= into the modified 



LLG equation leads to 0o : 
The total energy, ffs -t- Hst, 
then given by j29l | 



± -{ha^/2eS^XK_L)js < 1. 
of this static domain wall is 
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Since the energy iJuni of the uniform state is zero inde- 
pendently of js, we find that -Edw < £'uni when js > 
^^/2KK±X (Fig. 1(b)). This critical value is of the 
same order of magnitude as jg'' for K± ^ K. Thus, the 
uniformly magnetized state is a 'false vacuum', which 
collapses into a multi-domain state via first order phase 
transition under spin current above a critical value. 

The domain formation may be understood from the 
spin wave around uniform ferromagentism. In fact, the 
spin wave energy, ri™', given by Eq. indicates that 
the typical wavelength of the spin wave excitation (deter- 
mined by the minimum of the spin wave energy) is given 
by fc-i -A^(V1 + nKS)/{hvs) for smaU js- As the 
jsin~^(vs/tic)). Hence, current increases, the wavelength becomes shorter and 
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the critical current is given by the condition it becomes 
comparable to the wall thickness. 

The above estimate of the critical current density ob- 
tained for a planar wall would be applicable also to the 
case of nucleation in films such as in Ref . 0| , at least for 
an estimate of order of magnitude. 

Let us estimate the critical current numerically. The 
magnitude of the current needed for nucleation in a wire 
is estimated as js ~ 3 x 10^'^ A/m^ by using the materials 
parameters of Co: a = 2.3 A, S ^ 1, J = 3 x 10"''° W, 
K = 6.2 X 10-24 J, (A = ^/jJK = 76 A), under the 
assumption that K± is mostly due to the demagnetizing 
field, K± = M^a^/fiQ = 3 x 10^^^ J. Recent experi- 
ment on a Co layer with a current through a point con- 
tact with domain suggests a domain wall formation for a 
maximum current den sity close to the contact region of 
j ~ 5 X 10^2-13 A^/j^2 This value is in rough 

agreement with our estimate. Other recent experiments 
on metallic and semiconducting pillar structures indicat- 
ing inhomogeneous magnetization reversal may also be 
due to domain (or domain-like) structure created by spin 
current [Hill 113. 
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FIG. 1: Schematic phase diagram under spin current js in 
the absence of pinning potential, (a) K± < K. Above j", 
uniform ferromagnetism collapses into multi-domain structure 
in which domain walls are flowing due to spin current. The 
threshold, js for "depinning" from K±^ is below . (b) 
K± ^ %K. Energy of the single-wall state (-Edw) is com- 
pared with that of the uniformly magnetized state, -Buni = 0. 
Multi-domain state here remains at rest. In the gray region, 
(is > js'^'""), the domain wall starts to flow but is unstable, 
suggesting a new ground state. 



